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ABSTRACT 

Barycentric corrections made to the timing of Kepler observations, necessitated by 
variations in light arrival time at the satellite, break the regular time-sampling of the 
data - the time stamps are periodically modulated. A consequence is that Nyquist 
aliases are split into multiplets that can be identified by their shape. Real pulsation 
frequencies are distinguishable from these aliases and their frequencies are completely 
recoverable, even in the super-Nyquist regime, that is, when the sampling interval 
is longer than half the pulsation period. We provide an analytical derivation of the 
phenomenon, alongside demonstrations with simulated and real Kepler data for S Set, 
roAp, and sdBV stars. For Kepler data sets spanning more than one Kepler orbital 
period (372.5 d), there are no Nyquist ambiguities on the determination of pulsation 
frequencies, which are the fundamental data of asteroseismology. 

Key words: asteroseismology - methods: data analysis - stars: oscillations - stars: 
variables: 5 Scuti - stars: variables: general 
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1 INTRODUCTION 

It is well known that all light emitted or received is Doppler 
shifted by the motion of the emitting or receiving body. 
In astronomy the Doppler shift from the emitting body is 
used to deduce radial velocity; the Doppler shift of the re- 
ceiving body, the Earth, is removed from the observations 
by correcting the times of observation to the Solar System 
barycentre. All of this is standard practice. The Doppler 
shift from emitting bodies has traditionally been used to 
deduce orbital motion for binary stars and for exoplanets 
with spectroscopic radial velocities, and with O — C anal- 
yses in the cas es of binary pulsars and o rbiting pulsating 
stars. Recently. IShibahashi fc Kurtz! (|2012l ) have shown how 
this orbital Doppler shift can be seen directly in the Fourier 
transform of light variations in binary pulsating stars and 
used to derive the mass function without the need of spec- 
troscopic observations. Now, we show here, both analytically 
and in practice, that the Doppler shift at the receiving body 
- in this case the Kepler Space Telescope - allows the unique 
identification of periodic frequencies in the emitting body 
without any Nyquist ambiguity. 

The Nyquist frequency of equally-spaced data repre- 
sents an upper limit on a range of frequencies over which 
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a Fourier transform is unique. It is well-defined as /N y = 
l/2At, where At is the sampling interval between two con- 
secutive points; this is also known as the "cadence" of the 
data. In perfectly equally-spaced data there is an infinity of 
frequencies that fit the data equally well. External physical 
constraints are mandatory to select the frequency range that 
is appropriate for pulsating stars. Commonly, the cadence is 
selected so that the known plausible frequency range for the 
type of star lies in the range [0,/Ny], and the infinity of 
ambiguous higher frequency Nyquist aliases are ignored as 
unphysical. 

However, it is not possible to adjust the observing ca- 
dence for Kepler mission data to accommodate the study of 
pulsating stars in different frequency regimes. The Kepler 
mission's prime goal of detecting Earth-like planets in the 
habitable zone led to a design where Kepler data are avail- 
able in two cadences: long-cadence (LC) at 29.43 min, and 
short-cadence (SC) at 58.85 s, with corresponding Nyquist 
frequencies of 24.469d _1 and 734.07 d _1 , respectively. The 
number of short-cadence observing slots is limited to 512 
by telemetry constraints, and most of those observing slots 
are needed for higher time resolution studies of exoplanet 
transits. With this limitation most classical pulsators ob- 
served by Kepler have used LC, even though many such pul- 
sators - such as S Set, roAp, sdBV, pulsating white dwarf 
and ft Cep stars - have pulsation modes with frequencies 
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that can exceed the LC Nyquist frequency. iMurphvl (|2012T ) 
showed some of the difficulties encountered with Nyquist 
aliases when studying S Set stars with Kepler LC data. 

Now, we show in this paper how the correction of the 
times of observations of Kepler data to the Solar System 
barycentre generates a time-dependent Nyquist frequency, 
and that this completely removes all Nyquist alias ambigu- 
ities in the amplitude spectrum. For data sets longer than 
one Kepler orbit, true pulsation frequencies can always be 
distinguished from all Nyquist aliases. The observing ca- 
dence for Kepler data is not a barrier to asteroseismic mod- 
elling of pulsation frequencies because of Nyquist ambiguity 
- there is none. In practice, of course, hig her frequencie s 
have reduced amplitude with LC data (see lMurphy||2012T ). 
and that affects the signal-to-noise ratio, which decreases 
with increasing frequency. 

The Kepler spacecraft is in a 372. 5-d, heliocentric 
Earth-trailing orbit. Four times per orbit the satellite must 
perform a roll to keep its solar panels facing the Sun, so 
the data are divided into Quarters (denoted Qn) and a brief 
gap in observations occurs. In addition to the data downlink 
that takes place at each quarterly roll, two more downlinks 
occur at approximately 30-d intervals in the middle of each 
quarter, effectively dividing the data into three 'months' per 
quarter (denoted Qn.m for SC dat43). These events repre- 
sent the most frequent gaps in Kepler data which are other- 
wise nearly continuous at their micromagnitude precision. 
Other gaps due to occasional safe-mode events have oc- 
curred, and we refer the reader to the Data Characteristics 
Handbook for more information. However, the important 
property of the observations with respect to the gaps is that 
observations are always taken at regular intervals according 
to the spacecraft's clock, which can be expressed as some 
integer product with At from some (fixed) arbitrary start 
time, i.e. t n = to + nAt. In this regard, the Kepler data may 
be described as equally spaced. 

Timing onboard the spacecraft is in Julian Date, but 
the spacecraft is in orbit around the solar-system's barycen- 
tre and as such the arrival times of photons from the Kepler 
field are shifted seasonally by up to ±200 s from those of 
the Solar System barycentre. Although data are sampled 
regularly onboard, the time stamps are subsequently modi- 
fied when converted into Barycentric Julian Date (BJD) and 
hence become irregularly sampled: no longer can the obser- 
vation times be represented in the form t„ = to + nAt for 
integer values of n. The effect on super- Nyquist frequencies 
has been observed bv lBaran et"all l |2012h , and described as a 
'smearing' of the frequencies. In this paper we show how the 
satellite motion results in multiplets being generated out of 
the Nyquist aliases, and how this can be used to distinguish 
these aliases from real pulsation frequencies. 



1 LC data are not separated into different files around these 
months like the SC data are, but the gaps are still present. 

2 available at http: / /keplergo. arc. nasa.gov/Documcntation.shtml 



2 ANALYTICAL DERIVATION 
2.1 Discrete Fourier transform 

2.1.1 General description 

Let us consider first a general case of discrete data {x(t n )}, 
which can be regarded as the following function of any value 
oft; 



x N (t) =^2x(t)S(t-t n ), 



(1) 



where x(t) on the right-hand side is a continuous function 
x(t) and 8(t — t„) is Dirac's delta function. It is well known 
that the Fourier transform of Kjv(t) is expressed as a con- 
volution of the Fourier transform of x(t) and the sampling 
window spectrum Wn(u): 



N+l 



F n (uj) = 



1 



N + l 



XN(t) exp(iujt) dt 



(F * Wn)(ui) 



C2) 



where 



/oo 
x(t) exp(iwt) dt (3) 
- oo 

is the Fourier transform of the continuous function x(t) and 



W n (lo) 



N 



1 N 

— ^exp(iwt„). 



(4) 



Here, we have used the following expression of the delta 
function: 



/oo 
exp(— iuit) duo. 
-oo 



(5) 



That is, the Fourier transform of discretely sampled data 
is equal to the convolution of the Fourier transform of the 
original continuous function and the window function. If we 
assume x(t) — cos(u}ot + cf>), i.e. pulsation with a single mode 
of angular frequency too, regarding this as a model of a pul- 
sating star, then 

F(w) = i{J(w+wo) + <5(w-w )}. (6) 

Thus the Fourier transform of the sample series Fn {uj) con- 
sists of a superposition of the shape of the window spectrum, 
shifted by ±cjo- 



2.1.2 The case of a uniform cadence 

If the sampling is taken with a uniform cadence with a time 
interval At, then t n = to + nAt, and 



W n {uj) 



1 sin{(iV + l)wAt/2} 



N+l sin(ajAt/2) 

The window spectrum |W0v(w)| has a series of sharp, high 
peaks at ui — nws, where 

2vr 



At 



(8) 



and n denotes integers (see Fig.[T]). Then Fn(u) has sharp 
peaks at uj — nus + loo ■ These apparent multiple peaks are 
the Nyquist aliases, and they look identical except for their 
frequencies. 
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(coAt/2)/j[ 

Figure 1. The sampling window spectrum Wm{^) in the case of 
sampling with a constant time interval. The abscissa is uiAt/(2ir), 
where At is the sampling time interval. The red curve shows the 
case of N = 100 data points, while the green curve shows the 
case of N = 10. The window function shows conspicuous peaks 
at u) = riids, where o;g = 2ir/ At is the sampling angular frequency 
and n denotes integers. 



The peak in the range of [0,ct>s/2] corresponds to the 
true angular frequency as long as 2wo < (here, we take 
ujo > 0). Hence, for a given angular frequency ojo, if we select 
the sampling rate At < (27r/aio)/2, we can identify uniquely 
the true frequency from the Fourier transform of the sam- 
pled data. This condition of At < (27r/tJo)/2 is a sufficient 
condition for unique identification of the true frequency, not 
a necessary condition. 

Now let us consider a case that the sampling rate is 
given in advance. In this case, if the true angular frequency 
is lower than ojs /2, the lowest frequency peak of the Fourier 
transform of the sampled data gives the true angular fre- 
quency. Note that this is again a sufficient condition, not a 
necessary condition. The Nyquist frequency is this frequency 
limit, under which we can identify uniquely the true angular 
frequency: 

/Ny = o — - 2AT (9) 



2 2tt 



The reason why the true frequency is not uniquely iden- 
tified without external constraints in the case of sampling 
with a uniform cadence is the fact that all the Nyquist 
aliases look identical and are indistinguishable. Inversely, 
if we could distinguish the individual Nyquist aliases, we 
would be able to identify uniquely the true frequency even 
if cjo > ws/2. We will see in the case of Kepler data this is 
indeed possible. 

2.1.3 Sampling with a periodically modulated time interval 

In the case of Kepler, observations are taken at regular in- 
tervals according to the clock on-board the spacecraft. Data 
sampling is made at 



t n = t + nAt, 



(10) 



where n = 0, N, and At is a constant. The long-cadence 
(LC) data have At = 29.4 min. But, since Kepler is orbiting 
around the Solar System barycentre, the resulting annual 




equinox 



Figure 2. A schematic picture showing geometric relations 
among the Earth, the Sun, and the star. The Sun is moving 
around the Earth on the ecliptic plane. The geocentric ecliptic 
longitude is measured from the vernal equinox. The star's position 
is (A, 0), and the geocentric ecliptic longitude of the Sun is \o(t). 
The difference in the path length is shown with a thick segment. 
The path difference is obviously given as acos(A — Xq) cos /3. 



variation in the distance between the stars and the space- 
craft leads to modulation in the phase of observed stellar 
pulsation, and the time stamps are converted to barycentre 
arrival times. As a consequence, the time stamp interval of 
the Kepler data is periodically modulated. 

To make the modelling simple, we ignore in this section 
the difference between the barycentre and the heliocentre. 
Also, we approximate Kepler's orbit as a circle with a radius 
of a = 1 au. Then, the arrival time of the light from the star 
to the spacecraft delays from that to the heliocentre by 



5t(t) = — cos f3 cos{A0(t) — A}, 
c 



(11) 



which is often called the R0mer delay, in honour of Ole 
R0mer, who discovered observationally the fact that the 
speed of light is finite by comparing the observed timing of 
the eclipses of Jupiter's sa tellite Io with the estimate made 
at his time |Sterkenll2005Tl . Here, A and /3 are the ecliptic 
longitude and latitude of the star, A©(t) is the geocentric 
ecliptic longitude of the Sun, and c denotes the light speed 
(see Fig.[2j) . With the present assumption, 



x Q (t) = n(t-t ), 



(12) 



where f2 = 2-7r/365 radd - and to is the vernal equinox pas- 
sage time of the Sun. It is instructive to write down here, 
for later use, the light time with the equatorial coordinates, 

<5t(t) = — cos A© cos a cos 5+sin A© (cos e sin a cos 5+sine sin 5), 

(13) 

where a and S are the right ascension and the declination of 
the star and e is the obliquity of the ecliptic. 

By taking account of this light time effect, the Kepler 
data are stored as a function of the arrival time at the he- 
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Figure 3. The time correction for the case of A = 0°,/3 = 0° 
(red) and for the Kepler field-of- view (a = 300°, S = 50°) (green) 
as a function of the Sun's ecliptic longitude Aq. 



liocentr<jf]: 



where 



f 0?n = t„ + TCOs(£lt n — A), 



cos P 



(14) 
(15) 



gives the amplitude of modulation of the time stamps of 
the Kepler data. The second term on the right-hand side 
of equation (|14p is the heliocentric time correction. Fig. [3] 
shows this time correction (the R0mer delay) for the field of 
A = 0°,/? = 0° (red) and for the K epler field-of- view (green). 
The case of j3 = 0° gives the maximum correction; r = 500 s. 
In the case of the Kepler field, r ~ 190 s. 



2.2 The window spectrum of periodically 
modulated sampling 

Our aim is to carry out Fourier analysis of the Kepler data 
of pulsating stars. Since the Fourier transform of the sam- 
ple series Fn(oj) is a superposition of the window spec- 
trum, the problem is essentially what the window spectrum 

JV 

Wn(oj) = {N + l) 1 ^^ exp(iutQ :n ) looks like in the case of 

n=0 

periodically-modulated sampling. The problem is then how 
to treat the terms exp(iwrcosf2t) and exp(iujT sin fit). 



2.2.1 Bessel coefficients 

These terms can be expressed with a series expansion in 
terms of Bessel functions of the first kind with integer order. 



3 Strictly speaking, the data are stored as a function of the 
Barycentric Dynamical Time (TDB) referenced to the Solar Sys- 
tem barycentre. The difference between the time referenced to 
the heliocentre and that to the Solar System barycentre is < 4 s. 
The main cause of this difference is the acceleration of the Sun 
due primarily to Jupiter and Saturn. 



With the help of Jacobi- Anger expansions 

oo 

cos(±£cosvp) = J (£) + 2^(-l)"J 2 „(C)cos2n99 (16) 

71=1 

and 

oo 

sin(±£cosv?) = ±2^(-l) n J 2 „+i(Ocos(2n+ l)tp, (17) 

Tl = 

we easily obtain the following relations: 



cos(i9 ± £ cos if) = 
and 

sin($ ± £ cos ip) = 



J2 (-ir-jy2n(£)cos(tf + 2nv?) 

n— — oo 

T^2n+i(£)sin{tf + (2n+l)v?}j (ig) 
J2 (-l) m {./2n(£)sin(tf + 2^) 

n — — oo 

TJ'2n+i(Ocos{ 1 9 + (2n + 1M}. (19) 



Here Jn(£) denotes the Bessel function of the first kind of 
integer order n. 



2.2.2 Analytic expression of the window spectrum 

By applying these relations, after somewhat lengthy but 
straightforward manipulations, we eventually get 

(iV+l) 2 |W/jv(a;)| 2 

^ sin{(iV + l)(a; + fcn)At/2} 2 

~ fcf^ sin{(w + fcn)At/2} 

OO oo 1 

+ J2k(uT)J 2k ,(LUT) 
_ k= — oo fc' = — oo 

x cos I 2{k - k')VL^-^- - 2{k - k')X j 

OO OO ' 

+ ^+l(w)^2k'+l(wT) 
fc — — oo k' = — oo 

X cos | 2{k - k')Q^^ - 2(k - k')\ | 



OO oo 



fc — — oo fc ; — — OO 

x sin [ {2uj + [ 2(k + k') + l]Q} 
- {2(k + k') + 1} A] 

sin{(JV + l)(w + kQ)At/2} 



NAt 



sin{(w + kQ)At/2} 



(20) 



where means summation for k' except for when k' = k. 

The Bessel coefficients of different orders are orthogonal 
to each other (see Fig-[4j . Hence, among the terms on the 
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Figure 4. The Bessel coefficients J n (£) with n = 0, 8 for £ = 
[0 : 5.5]. The vertical lines show £ = nuigr (n = 0, 8), at which 
the window spectrum has sharp multiplct peaks (see text). 



Table 1. The window spectrum. 
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uj s ± 2U 




2 





2ui s 


Jo(2u;st) 




±1 


2uj s ± O 


■/i(2w s t) 




±2 


2oj s ± 20 


J 2 (2l> s t) 




±3 


2ui s ± 3f2 


M2u, s t) 


3 





3cj s 


Jo(3ust) 




±1 


3o;g i 


J 1 (3ui s t) 




±2 


3w s ± 2Q 


J 2 (3ujst) 




±3 


3o)s ± 3f2 


J3(3wst) 




±1 


3u s ± 40 


J4,(3u>st) 



right-hand side, the first term is dominant. Then, 
1 



N+l 



x ^2 m^) 



sm{(N + l)(u + kQ)At/2} 



sin{(o; + fcft)Ai/2} 



(21) 



This means that the window spectrum | Wn(oj) | consists of a 
singlet sharp peak at uj = and multiplets of sharp peaks at 
uj = nuis + kQ, where n and k are integers and ujs = 2n/ At. 
The amplitude of each peak is |Jfc(nwsr)|, since <C 
(see Tabled]). 

It should be stressed here that |Wjv(a;)| is highly de- 
pendent on frequency uj through £ = lot, the argument of 
the Bessel function. For f < 1, the dominant term is only 
Jo(0 ~ 1— £ 2 /4 and all the higher-order Bessel functions are 
negligibly small. However, with the increase of £, the first- 



order Bessel function Ji(£) comes larger as Ji(f) ~ f/2, 
and eventually it becomes larger than Jo(£) around £ ~ 1.5. 
But the higher-order terms Jfc(£) for > 4 are still neg- 
ligibly small. This means that |Wjv(w)| looks like a triplet 
with an equal spacing of Q if uj < (2t) _1 . With the further 
increase of f, the second-order Bessel function J2(£) also 
becomes non-negligible. As a consequence, |VFjv(t<j) comes 
to look like a quintuplet. This tendency continues further: 
the high-order Bessel functions become non-negligible with 
the increase of uj. Therefore, the individual Nyquist aliases 
are distinguishable from each other. This characteristic is 
different from the case of sampling with a uniform cadence. 



2.2.3 Nyquist frequency is no longer the detection limit 



We have shown in Sec. 12. 1.11 that the Fourier transform of 
the discretely sampled data a;jv(t) is given as the convo- 
lution of the Fourier transform of the original continuous 
data x(t) with that of the window spectrum Wn(uj). Hence 
the Fourier transform of discrete data with a periodically- 
modulated sampling interval essentially consists of multi- 
plets around the Nyquist aliases of Wn{uj) at nujs ± ujq, 
where n denotes integers. It should be noted here that the 
shape of the multiplets is independent of the frequency of 
the mode loq, as long as the amplitude and the frequency 
of the mode are stable. It should also be stressed that only 
the peak at to — loq with n — is a single peak, while the 
other Nyquist aliases are multiplets; nuos ± loo ± kQ, where 
n 7^ and k denotes integers (see Table[2]). Hence, by finding 
the unique singlet, we can distinguish the true pulsation fre- 
quency, uo, from its Nyquist aliases. Note also that this fea- 
ture is true not only for the case of uj < los /2 but also for the 
case of uj > ujs/2. That is, a pulsation frequency higher than 
the Nyquist frequency, which is still defined as l/(27rAt), is 
uniquely determined when the sampling rate is periodically 
modulated, irrespective of the number of multiples of the 
Nyquist frequencies crossedQ The Nyquist frequency is no 
longer the upper limit of frequency determination in such a 
case. This has been mathematically proved here. 

2.3 Window spectrum for the Kepler LC data 

In the case of the Kepler field, r ~ 190 s ~ 2.20 x 10 -3 d. 
For the range of w/(2ir) = [0,200] d _1 (uj = [0,1.26 x 
10 3 ] radd -1 ), ujt = [0, 2.77] rad. Fig. 0] shows the Bessel 
coefficients J„(lot) with n = 0, 8 for ujt — [0, 5.5] . From 
this figure, it is obvious that J n (uJT) with n ^ 8 are negligi- 
bly small in this range. 

Let us consider the case of LC of Kepler. If we take 
At = 30min = 2.1 x 10~ 2 d, then uj s = 2n / At = 3.016 x 
10 2 radd _1 . Then lust = 0.663 rad. As for the modulation, 
Q = 2tt/365 rad d" 1 = 1.72 x 10~ 2 radd~ 1 , then Q/uj s = 
5.69 x 10~ 5 . Parameters for the Kepler LC sampling are 
given in Table[3] 



4 For truly irregularly (randomly) sampled data, the lowest fre- 
quency above which real frequencies are truly indistinguishable 
from aliases is eq ual to the reciprocal of the accuracy to which 
time is measured For Kepler this is around 10 d , 

so we can distinguish aliases from real frequencies for all frequen- 
cies within the star that are physically meaningful. 
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Figure 5. The dependence of (jV + Jki^r) sin{(N + l)(w + fcf2)At/2}/ sin{(o; + fcfi)Ai/2} | on w/w s > for fc = 0, ...,8, where 

Log = At/(2n). The color codes, corresponding to the order of the Bcsscl coefficients k, are in common with Fig. [4] The parameters are 
At = 30min; ui s = 3.016 X 10 2 radd" 1 . O = 2tt/365 rad d" 1 ; ft/i^s = 5.69 X 10" 5 . r = 190s = 2.20 X 10~ 3 d; oj s t = 0.663 rad; AT = 10 5 . 



Fig.[S]shows the dependence of | Jk(^ t) sin{(A r +l)(w + 
kQ)At/2}/ sin{(u; + kQ)At/2} \ on oj/uj s for fc = 0, ...,8. 
Here, we assumed N = 10 5 , that is, the observational time 
span was assumed to be ~ 5.75 yr. In such a case, the high 
peaks are so sharp that their self-cross terms dominate over 
the other cross terms. As seen in Fig.O the high peaks ap- 
pear as multiplets around lj = mos, where ljs = 2n/At and 
n denotes integers. The multiplets have the equal splitting 
of the orbital frequency, SI (f2/ws = 5.69 x 1CT 5 in Fig. [5]). 

Strictly speaking, each multiplet has infinite peaks. 
However, the tiny amplitude peaks are practically unde- 
tectable, and each multiplet looks like a finite number of 
components; the first Nyquist alias looks like a triplet, the 
second looks like a quintuplet, and so on. The degree of the 
apparent multiplicity becomes higher with the increase of n, 
because the Bessel coefficients with n/0 become important 
with the increase of oj. 

The amplitudes of the peaks in the individual multiplet 
are determined by ojst. Since r ~ 190 s ~ 2.20 x 10 -3 d, 
ojst = 0.663, 2uj s t = 1.326, 3ui s r = 1.989, 4w s r = 2.652, 
5uj s t = 3.315, 6uj s t = 3.978, 7cj s t = 4.641, and 8lu s t = 
5.304. The Bessel coefficients at these values are graphically 
seen in Fig.(4] 



2.4 Multiplets of the Kepler LC data 

It should be stressed here again that only the peak at uj = loq 
with n = is a single peak. The aliases with frequencies 
±ujs ±<^o look like a triplet, while those with 2cjs ±<^o are a 
quintuplet. The aliases associated with 3ws, 4ws, 5o;s, and 
6ojs look like a septuplet, a nonuplet, a undecuplet, and a 
tredecuplet, respectively. 



3 EXAMPLES WITH SIMULATED DATA 

3.1 Importance of a large observational time span 

With the real Kepler data we are limited in our observa- 
tional time span by the length of the mission thus far, as- 
suming that the star we wish to study has been observed 
continuously since the mission began - 3 complete orbits of 
the satellite have elapsed since then at the time of writing. 
However, with simulated data we are not so restricted. In 
Fig.[S]we show the importance of a large observational time 
span on resolving the multiplets into which Nyquist aliases 
are split. Although these simulated data are without noise, 
it is clear that when one orbit (four mission quarters) of 
data is available, aliases are no longer single peaks, but are 
split into sets of partially-resolved peaks and can thus be 
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Table 2. Fourier transform of the data sampled with a periodi- 
cally modulated interval. 
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Table 3. Parameters for the Kepler LC sampling. 
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distinguished from real peaks (which are neither split nor 
similarly distorted). With two or more orbits, an obvious 
multiplet emerges, but the full-width-half-maximum of each 
peak is still large compared to the separation of the multiplet 
peaks, hence the window pattern of each peak can interfere 
with its neighbours. Increasing the number of orbits further 
sharpens the peaks and resolves the pattern for comparison 
with the analytical patterns we derived, shown in Fig. [5] 



3.2 Effect of varying co-ordinates 

It is the barycentric time corrections that cause the split- 
ting of Nyquist aliases. For a satellite that orbits around the 
Solar System barycentre in the ecliptic plane, these time 
corrections are determined by the ecliptic latitude of the 
target star, because this determines the difference in light 
arrival time at the satellite and the barycentre. We show how 
the multiplet shape changes with varying ecliptic latitude in 
Fig. [7] In the more familiar celestial co-ordinates, both right 
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Figure 6. Simulated data showing the frequency resolution 
achieved with different numbers of orbits of the satellite. For short 
observational time spans, the window pattern of each member of 
the multiplet distorts the shape of its neighbours. Calculated for 
Kepler field co-ordinates, at the first Nyquist alias of a 10 d — 1 
pulsation frequency. 



© 0000 RAS, MNRAS 000, 000-000 



8 S.J. Murphy et al. 



ascension and declination affect the ecliptic latitude, so we 
expect the arrival times and thus the multiplet shape to 
change when either co-ordinate is varied. Since neither of 
these coordinates changes greatly over the Kepler field, the 
multiplet shape is almost identical, too. 



4 SUPER-NYQUIST ASTEROSEISMOLOGY 

4.1 Demonstration with real data 

We now compare our theory with real Kepler data by ex- 
amining the case of a high-amplitude 5 Set (HADS) star. 
These stars are the main-sequence counterparts to the clas- 
sical Cepheids, having large amplitudes and normally few 
pulsation modes - often the fundamental radial mode alone 
or with the first radial overtone. Sometimes harmonics of 
the frequencies of these modes are seen, as with the case 
we present here: that of KIC 6861400. We used the Ql-9 
LC data, which have been processed with the PDC-MAP 
pipeline jStumpe et alj|2012l : ISmith et al.l 120121 ). The star 
features three independent mode frequencies and many har- 
monics of the two mode frequencies with the highest ampli- 
tudes. We removed ('pre- whitened') all the harmonics from 
the data, and kept only the highest-amplitude mode of the 
three, leaving a single, high-amplitude peak to make our 
example. No data points were deleted. 

As expected, there is one alias in each frequency range 
[n/Ny, {n+ 1)/Ny] (upper panel, Fig.|HJ), for integer n. Even 
values of n, counting from n = at 0d~ , are multiples of 
the sampling frequency. In Fig. [HI two peaks straddle each 
of these, and are identical in the Fourier transform. For this 
reason, only the peak to the left of each multiple of the 
sampling frequency is presented in the lower panels for a 
closer look. The multiplet structure is clearly present in all 
but the real peak, which is represented by a unique singlet 
peak. The spacing in the multiplets is approximately equal 
to the orbital frequency of the satellite, but due to the small 
number of orbits covered (2.25), the frequency resolution 
is insufficient to separate fully the window patterns of the 
members of the multiplets. The window pattern of each peak 
of the multiplet thus interferes with neighbouring multiplet 
members, and distorts the shape somewhat, as we showed 
in Fig.(S] 

As the frequency increases when we study consecutive 
aliases of higher n, the argument to the Bessel functions 
is also higher and the higher-order (higher k) terms be- 
come more significant. As such, the number of peaks seen 
increases. The amplitudes of the more distant members of 
the multiplet increase - we begin to see the expected quin- 
tuplets, septuplets and so on (cf. § 12. 4p . as these peaks' am- 
plitudes climb above the noise level. The amplitude ratios of 
the peaks with respect to each other also change. In panel 
a, the central peak is dominant, but it is absent in panel c, 
where only the window pattern of the first order sidelobes 
can be seen in its place. The amplitudes of the Bessel or- 
ders that determine these peaks' amplitudes were shown as 
a function of frequency in Fig.[5] 

The frequency dependence of the shape of the multiplet 
is on the sampling frequency straddled, only. The multiplet 
generated at the first Nyquist alias of a pulsation frequency 
at 20 d _1 looks exactly the same that of a 10-d _1 frequency, 
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Figure 7. Variation in multiplet shape for changing ecliptic lat- 
itude. Note that, because the ecliptic is inclined to the equator, 
varying cither right ascension or declination will change the eclip- 
tic latitude. At an ecliptic latitude of 90°, a single peak is seen 
because there is no longer any difference in light arrival time. 
Calculated for 20 orbits, at the first Nyquist alias of a 10 d 
pulsation frequency. 
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Figure 8. The upper panel shows the frequency distribution of Nyquist aliases. In the lower panels, which are all to the same scale in 
frequency and amplitude, the left column shows the real peak, as labelled in the upper panel, and its window function. Panels a-d show 
aliases (also labelled in the upper panel), and their relationship to the sampling frequency and the real frequency. The pattern of the 
number of peaks, p, as a function of the coefficient of the Nyquist frequency straddled, n, is apparent: p = In + 1. As such, there are 
three peaks at / = /g + /r, three at /g — /r, five each at 2/g ± /r and so on, as shown in Fig. [5] The example shown is KIC 6861400, 
using Ql-9 PDC-MAP data. 



and so on. Only the frequency of the alias changes, not the 
nmltiplet shape. 

4.2 Distinguishing real peaks from Nyquist aliases 

The difficulty in observing S Set stars with Kepler is that the 
typical pulsation frequency range of S Set stars spans both 
sides of the LC Nyquist frequency. Those pulsation frequen- 
cies above the LC Nyquist frequency are aliased to lower fre- 
quencies in the Fourier transform, meaning that they over- 
lap with the real peaks. Distinguishing these aliases from 



the real peaks is impossible with data that are truly equally 
spaced in ti me (cf. § 12.1.20 . even if gaps exist in the data 
(|Koenll2010r ). It has been argued that SC data are required 
to study the S Set stars so that this problem of overlap- 
ping real and aliased frequencies is removed, because the 
SC Nyquist frequency is so much higher than the LC one. 
With a periodically-modulated sampling, however, one can 
use the shape of a peak to determine whether it is real or 
aliased. Fig.[9]shows how a real peak exists as a single peak in 
the Fourier transform, but a Nyquist alias will be a multiplet 
(resolved or not) split by the satellite's orbital frequency. 
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Figure 9. In the left panel, the SC (top) and LC (bottom) data are displayed. LC Nyquist frequency and sampling frequency are 
represented with vertical orange lines. The LC data lying above their Nyquist frequency are coloured magenta; those below are red. 
There arc Nyquist aliases in the LC data where real peaks have been aliased across the Nyquist frequency. By comparison with the 
SC data, we can sec these aliases have no real partner in the SC data. An example of this is the alias peak labelled 'a', corresponding 
to the real peak labelled V — there is no peak in the SC data at the same frequency as peak 'a'. We zoom in on 'a' and 'r' from the 
LC data in the right panel. The real peak takes on the typical sinc-function shape as expected in the Fourier transform. The effect of 
periodically-modulated sampling, however, is that the aliased peak is split into a multiplct, provided the time span of the data is at least 
one complete Kepler orbit. The example given is the S Set star KIC 8590553. SC data cover Q4.3 only, but LC data are from Q0— 9. 



Use of this technique to distinguish Nyquist aliases 
opens up asteroseismic possibilities for many Kepler tar- 
gets that have been observed in LC only. With data sets 
that have time spans greater than one Kepler orbital period 
there is no Nyquist ambiguity in selecting the true pulsation 
frequencies for all types of pulsatin g stars. While there are 
still benefits of SC over LC data |Murphv||20T3) . the SC 
data are not needed to resolve Nyquist aliases. 



4.3 Beyond the Nyquist frequency 

The rapidly-oscillating Ap stars are a population of chem- 
ically peculiar A-type stars that pulsate with frequencies 
much higher than S Set stars - frequencies between about 
70 and 250 d _1 . As such, when viewed with LC data, these 
pulsations are many multiples of the Nyquist frequency away 
from falling into the range [0, fyy]- It is still possible, how- 
ever, to identify the true frequency (/r = ujq /2tt) of the pul- 
sations with LC data, even though the sampling period is 
much longer than the pulsation period, if the Fourier trans- 
form is calculated over a frequency range that includes /r. 

We present an example in Fig. [TO] using the star 
KIC 10195926, whose roAp pulsations have periods among 
the longest known and whose frequency spectrum features 
an oblique dipol e mode at ~ 84 d~ that is rotationally split 
into a septuplet (|Kurtz et al.l201lT ). We have run a high-pass 
filter on the data to prewhiten the low frequency content. 
Thus the Ql-11 LC data presented contain only the high 
frequency pulsation. With this figure we present a few final 
examples of the application of the technique to LC data. 

Firstly, it is clear that despite being substantially higher 
in frequency than the Nyquist frequency, the real peak is 
recoverable and frequency analysis can be conducted with- 
out SC data. Secondly, peaks with the same relationship to 
the sampling frequency, i.e. the pairs matching nfs ± /r, 
are split into multiplets of the same shape, e.g. the multi- 



plets corresponding to 2/s — /r and 2/s + fn have the same 
shape, as is seen by comparing the lower panels a and /. 
Thirdly, even though the peaks have low amplitudes (being 
only about an order of magnitude higher in amplitude than 
the noise), it is a testament to the quality of Kepler data 
and confirmation of the validity of our theory that we can 
see so many orders of each multiplet. Specifically, panel e of 
Fig. [10] corresponds to the fifth Nyquist alias. We thus ex- 
pect 5n + 1 = 11 equally-spaced peaks to be visible. We plot 
blue circles at intervals of (1/372.5) d _1 , i.e. the orbital fre- 
quency, from the central component of the multiplet. Each 
multiplet member can be seen to lie at its predicted position, 
beneath a blue circle. There are 11 blue circles, and all but 
the left-most one corresponds to a clear multiplet member. 
This left-most member is of low amplitude, barely above the 
noise level and is therefore not resolved from another noise or 
window-pattern peak. Its counterpart, under the right-most 
blue circle, is identifiable. If we had a longer observational 
time span the peaks would sharpen and cross-talk between 
their window patterns would decrease, likely leading to the 
final unresolved member becoming more distinguished. 



4.4 Distinction from other modulation 

Periodic amplitude and frequency modulation of a stellar 
pulsation signal both generate multiplets in the frequency 
domain. The critical factor that distinguishes the Nyquist 
aliases we describe here from frequenc y modulation multi- 
plets in binary stars as described by Shibahashi fc Kurt j 
(|2012f ) . or amplitude modulation multiplets that are seen 
in, for example, roAp stars and Blazhko RRLyrae stars, is 
that the Nyquist alias multiplets are split by exactly the 
Kepler orbital frequency (1/372.5 d), within the frequency 
resolution. It is exceedingly unlikely that other targets will 
have the same modulation frequencies, hence show the same 
splitting. In the unlikely event that such a coincidental fre- 
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Figure 10. In this example of the roAp star KIC 10195926, a real multiplct of peaks appears around 84 d~ 1 . and is shown zoomed in the 
top-right panel. The many aliases thereof are visible along with the real multiplet in the top left. The Nyquist frequency and its integer 
multiples are indicated as vertical orange lines. In the lower panels, we demonstrate the effect of periodically-modulated sampling on the 
shape of the highest-amplitude peak in each dipole-mode multiplet, all to the same scale. The real peak (far-left, top) has the typical 
sinc-function shape, even in these LC data and well beyond the Nyquist frequency. It can be compared to the window function in the 
far-left, bottom panel. The remaining lower panels, a-f, correspond to labels in the upper panel, and are in frequency order, clockwise. 
The data are Ql-11 LC data. See the text for further explanation. 



quency were found in some star, then careful examination of 
the sidelobe amplitudes and phases would clearly distinguish 
the cases. Section 2 provides all the necessary information, 
should amplitude and phase need to be examined in this 
rare case. 



5 APPLICATION TO THE KEPLER SC DATA 

In the case of Kepler SC data, the sampling time interval 
is much shorter (at 58.9 s) than in LC mode. Hence in SC 



we have: ljs = 9.217 x 10 3 radd , and then ljst = 20.3 rad 
and Q/u s = 1.866 x 10" 6 . The SC Nyquist frequency is 
734.07 d- 1 (=8496.18 /iHz). 

Since the argument of the Bessel coefficients, ljst, is 30 
times larger in the case of Kepler SC than that of LC, the 
multiplet becomes much higher order even in the case of the 
lowest Nyquist alias, and also the amplitude of each peak 
shrinks as the power is distributed into more peaks. In the 
case of the lowest Nyquist alias, at ui — cos , the multiplicity 
is as high as 50, and the amplitude is reduced to about 
20 per cent of that of the true peak since the amplitude of 
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Figure 11. The window spectrum for the lowest Nyquist alias at 
Ul = of the Kepler SC data. 



the Bessel coefficients Jfc(£) asymptotically decreases as ~ 
[2/(7r^)] 1//2 . This means that, in the case of Kepler SC, each 
of the Nyquist aliases in the power spectrum look like a 
'forest' of peaks, while the true peak is a singlet five times 
higher in amplitude than the 'forest'. Hence the singlet true 
peak looks obviously and conspicuously different from the 
aliases, irrespective of whether the true frequency is higher 
or lower than the Nyquist frequency. The true peak should 
be more easily distinguished than in the LC case. 

Fig.[TT] shows the window spectrum for the lowest 
Nyquist alias at uj = uis in the case of N= 10 5 . Most of the 
peaks therein are as short as 20 per cent of the true peak, 
whose amplitude is normalised as unity, or much shorter. 
The multiplet would look like an unresolved, broad-band 
plateau if the resolution were lower. 

As a real example, we show in Fig.[T5] the alias peak 
in the SC data of our previous roAp star example. The up- 
per panel demonstrates the appearance of the oblique dipole 
mode in the frequency range beyond the Nyquist frequency. 
One can see that the Nyquist aliases are at least six in num- 
ber, where the seventh is at the noise level. Each peak looks 
like a thick pillar. That these are aliases is obvious because 
there is a forest where there would be a single peak if we 
were looking at the real pulsation frequency. The dashed red 
lines in the upper panel show the region that is plotted in 
the lower panel. 

As seen in this zoomed-in plot, the shape of the 'forest' 
matches reasonably well with the expected window spectrum 
shown in Fig. 111! in particular, for the right half (on the 
left some peaks are unresolved in the real data) the relative 
amplitudes of each peak fit well with the expectation, as 
does the amplitude ratio of the alias to the true peak: the 
multiplet has amplitudes that range from around 30/imag 
for the unresolved part, to 45 /xmag for the right hand side, 
while the amplitude of the true peak in the Q6-12 SC data 
is 168 /imag. The anticipated amplitude reduction to around 
20 per cent is confirmed. 

We provide one more example: that of aliasing of a pul- 
sation frequency above the SC Nyquist frequency. We use 
the sub-dwa r f B s tar KIC 10139564 that was examined by 
iBaran et af] (|2012h . In this star, non-linear combinations of 
mode frequencies lead to combination frequencies above the 
SC Nyquist frequency. It should be stressed that these fre- 
quencies are still r eal frequencie s that describe the light vari- 
ations of the star. IBaran et al.l correctly identified aliases of 
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Figure 12. The Nyquist alias of Kepler SC data of the roAp star 
KIC 10195926. The frequency of this alias is / s - / R . The upper 
panel shows what the mode looks like as a whole. The region 
delimited by the dashed red lines is shown in the lower panel, 
zooming in on the 'forest'. 



a few of these in their figure 20. The aliases identified were of 
the form /r — /s, and their Fourier calculations were based 
on 626 000 data points covering 462.5 d of SC data. Here, we 
extend the analysis to 10 6 data points covering 739 d. 

In addi tion to the combinat ion frequencies at 995.43 d _1 
reported bv lBaran et alJ l)2012h . lies a Nyquist alias centred 
at 995.31 d" 1 . This frequency to which this Nyquist alias be- 
longs was pre-whitened by those authors (identified there as 
/is = 472.855 d _1 ) and so its Nyquist alias does not appear 
in their figure 20. We presen t the situatio n in Fi g. [l3l 

Panel a shows '/is' of IBaran et all (|2012l ). In lower 
panels only, fig has been pre-whitened. Its Nyquist alias, 
fs — /is, is the broad forest in panel b, of substantially 
lower amplitude than the real peak, and in accordance with 
our theoretical expectations regarding both location in fre- 
quency and reduced amplitude. Also visible on the right of 
panel b are some low-amplitude peaks that are combination 
frequencies in this sdB star. If fig, is prewhitened, its alias 
multiplet disappears from panel b; what remains is visible in 
panel d and includes those combination frequencies. Finally, 
those combination frequencies have Nyquist aliases of their 
own, but since those (real) combination frequencies are so 
low in intrinsic amplitude, their Nyq uist aliases ju st con- 
tribute to the noise left in panel c - IBaran et al.l noticed 
and commented upon this increased noise, describing it as 
'low-amplitude complex signal'. 

What we are seeing is two sets of frequencies (/is and 
the combination frequencies), that lie almost equidistant 
from and on opposite sides of the SC Nyquist frequency, 
such that the aliases of one set fall right next to the real 
frequency of the other set, and vice versa. Panel 6 in par- 
ticular highlights how easily these Nyquist aliases can be 
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Figure 13. Horizontally adjacent panels have the same amplitude scale; verticall y adjacent panels co ver the same frequency region. All 
panels are to the same scale in frequency. The tall peak of panel a is the '/i 8 ' of iBaran etahl j2012h . and is pre-whitened in the lower 
panels, only. The Fourier transform is calculated on Kepler Q5-12 SC data, processed with pipeline version 8.0/8.1. Further explanation 
is provided in-text. 



distinguished from real peaks, even (and indeed especially) 
in this very high frequency regime. 



6 CONCLUSIONS 

The barycentric time corrections applied to Kepler data 
break the regularity of the time interval between consecutive 
observations - Kepler data are not equally spaced in time. 
For data spanning at least one orbital period, periodically- 
modulated sampling causes multiplets to be generated out 
of Nyquist aliases, whereas real peaks remain as singlets, ir- 
respective of whether they lie above or below the Nyquist 
frequency. Multiplicity of the Nyquist aliases, along with 
relative amplitudes between multiplet members, depends on 
the number of multiples of the Nyquist frequency crossed, 
the amplitude of the sampling modulation, and the obser- 
vational time span of the data. 

We have theoretically derived the expected shapes of 
the alias multiplets and shown that the observed patterns 
are in agreement with the theory. We investigated the shape 
as a function of the number of multiples of the Nyquist fre- 
quency crossed, and obtained the following results: In real 
Kepler LC data, where noise is small but none the less 
present, the Nyquist aliases are split into multiplets whose 
members number 2n + 1, where n is the number of multiples 
of the Nyquist frequency crossed. In the SC case, multiplicity 



is substantially higher and amplitudes are further reduced. 
Multiplet members are equally separated in frequency by 
the orbital frequency of the satellite in both cadences. The 
shape of the multiplet is independent of the actual pulsation 
frequency. 



The resulting amplitude reduction upon aliasing pre- 
cludes this method from being applicable to frequencies of 
very low inherent amplitudes, such as those of solar-like os- 
cillators - the amplitude reduction leaves no significant sig- 
nal. The major consequence of our findings is that LC data 
may be used to study <5 Set stars and roAp stars with real 
pulsation frequencies above the LC Nyquist frequency. This 
is also true using SC data for much higher frequency pul- 
sators such as sdBV stars and pulsating white dwarf stars 
that have frequencies above the SC Nyquist frequency. The 
requirement is that continuous data sets are needed that 
have time spans greater than the orbital period of Kepler. 
This is the case for many hundreds of stars observed by 
Kepler in LC that previously were thought to have insur- 
mountable ambiguities in their frequencies. 
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